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It may be possible to construct a laser interferometer gravitational wave antenna in space with 
hrms ~ 10~^^ at f ^ O.lHz in this century. We show possible specification of this antenna which we 
call DECIGO. Using this antenna we show that 1) typically 10^ (10* ~ 10®) chirp signals of coalescing 
binary neutron stars per year may be detected with S/N ~ 10*. 2) We can directly measure the 
acceleration of the universe by ten years observation of binary neutron stars. 3) The stochastic 
■ gravitational waves of Qaw ^ lO"'^" predicted by the inflation may be detected by correlation 

' analysis for which effects of the recent cosmic acceleration would become highly important. Our 

, formula for phase shift due to accelerating motion might be also applied for binary sources of LISA. 
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I. INTRODUCTION 



There are at least four methods to detect gravitational waves. They are ;1) Resonant type antenna covering ~kHz 
band; 2-a) Laser interferometers on the ground covering 10Hz~ kHz band; 2-b) Laser interferometers in space Hke 
CsJ ■ LISA covering 10^*^ ~ lO^^Hz band;. 3) Residuals of pulsar timing covering ~ lO^^Hz band; 4) Doppler tracking 
of the spacecraft covering 10^"' ~ lO^^Hz band. It is quite interesting to note that little has been argued on possible 
detectors in 10~^ ~ lOHz band. In this Letter we argue in §2 possible specification of such a detector which we call 
' DECIGO (DECi hertz Interferometer Gravitational wave Observatory). In §3 we argue that the direct measurement 
QQ of the acceleration of the universe is possible using DECIGO. §4 will be devoted to discussions. 

o 

II. SPECIFICATION OF DECIGO 

The sensitivity of a space antenna with an arm length of 1/10 of LISA and yet the same assumption of the technology 
Q ■ level, such as a laser power of 1 W, the optics of 30 cm, etc. will be 4 x 10~^ ^Hz~''^^^ around 0.1 Hz in terms of strain, a 
H factor of 10 better than the planned LISA sensitivity around 0.1 Hz (see also http: / /www. physics. montana.cdu/maggie 



^ for a project named MAGGIE around this band). The sensitivity could be improved by a factor of 1000 for the next 
. . [ generation of a space antenna with more sophisticated technologies such as implementation of higher-power lasers and 
^ ■ larger optics in order to increase the effective laser power available on the detectors, and thus to reduce the shot noise. 

The ultimate sensitivity of a space antenna in the far future could be, however, 3 x 10~^^ around 0.1 Hz in terms 
5_j of strain, assuming the quantum limit sensitivity for a 100 kg mass and an arm length of 1/10 of LISA. We name 
Ci this detector DECIGO. This requires an enormous amount of effective laser power, and also requires that the other 
noise sources, such as gravity gradient noise, thermal noise, practical noise, etc. should be all suppressed below the 
quantum noise. Here we assume that such an antenna may be available by the end of this century, although we note 
that within the next five years or so NASA will begin serious discussions of a follow-on to the planned NASA/ESA 
LISA mission, so DECIGO technology may be achieved sooner. Note here that when the pioneering efforts to detect 
the gravitational waves started in the last century using resonant type detectors as well as laser interferometers, few 
people expected the present achievement in resonant type detectors such as IGEC(bar) ||] and in laser inteferometers 
such as TAMA300 ||, LIGO, GEO600, and VIRGO (for these detectors see [g). Therefore all the experimentahsts 
and the theorists on gravitational waves should not be restricted to the present levels of the detectors. Our point 
of view in this Letter is believing the proverb " Necessity is the mother of the invention" so that we argue why a 
detector like DECIGO is necessary to measure some important parameters in cosmology. 

The sensitivity of DECIGO, which is optimized at 0.1 Hz, is assumed to be limited only by radiation pressure noise 
below 0.1 Hz and shot noise above 0.1 Hz. The contributions of the two noise sources are equal to each other at O.lHz, 
giving the quantum limit sensitivity at this frequency. The radiation pressure noise has a frequency dependence of 
(X /^^ (in units of Hz^^^^) because of the inertia of the mass, while the shot noise has a dependence of approximately 
cx (in units of Hz^^^^) because of the signal canceling effect due to the long arm length. In figure 1 we show 
sensitivity of various detectors and characteristic amplitude he for a chirping NS-NS binary at z = 1. 
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III. DIRECT MEASUREMENT OF THE ACCELERATION OF THE UNIVERSE 



Recent distance measurements for high-redshift supernovae suggest that the expansion of our universe is accelerating 
I which means that the equation of the state of the universe is dominated by "dark energy" with p + 3p < 0. 



SuperNova / Acceleration Probe (SNAP, http://lbl.gov ) project will observe ~ 2000 Type la supernovae per year up 
to the redshift z ~ 1.7 so that we may get the accurate luminosity distance cIl^z) in near future. Gravitational wave 
would be also a powerful tool to determine dL{z) 

From accurate dL{z) one may think that it is possible to determine the energy density p{z) and the pressure p{z) 
as functions of the redshift. However as shown by Weinberg and Nakamura & Chiba 1^, p{z) and p{z) can not be 
determined uniquely from d^lz) but they depend on one free parameter Qko (the spatial curvature). 

Recent measurement of the first peak of the anisotropy of CMB is consistent with a flat universe {ilko = 0) 
for primordially scale-invariant spectrum predicted by slow-roll inflation j^] under the assumption of A cosmology. 
However it is important to determine the curvature of the universe irrespective of the theoretical assumption on 
the equation of the state and the primordial spectra also. In other words an independent determination of Qko is 
indispensable since is by far the important parameter. As discussed in the direct measurement of the cosmic 
acceleration pO| ] can be used for this purpose. Here we point out that the gravitational waves from the coalescing 
binary neutron stars at z 1 observed by DECIGO may be used to determine flko- Even in the worst case the 
redundancy is important to conflrm such an important flnding as the dark energy. 

A. Cosmic Acceleration 

We consider the propagation of gravitational wave in our isotropic and homogeneous universe. The metric is given 
by ds^ — —dt^ A- a{t)'^{dx^ + r{x)^{d9^ + sin^ 9d(f>^)), where a{t) is the scale factor and a{t)r{x) represents the angular 
distance. The relation between the observed time of the gravitational waves to at x = and the emitted time t,, at 
the fixed comoving coordinate x is given by ° = x = const. Then we have dt^/dt^ = ao/a^ = (1 -f z) and 

dHo 



dtl 



(1 + z)a-\dta{to) - dtaite)) = gcos{z) = (1-t- z)((l -f z)Ho - H{z)), (1) 



where H[z) is the Hubble parameter at the redshift z and Hq is the present Hubble parameter. For an emitter at the 
cosmological distance z > 1 we have gcos{z) ~ 0{t^^) where t^ is the age of universe to ~ 3 x 10"'^'^ sec. From above 
equations we have At^ — Ate(l-I- z) + l^t^^ -I- • • • , where At^ and Atg are the arrival time at the observer and the 

time at the emitter, respectively. When we observe the gravitational waves from the cosmological distance, we have 

Mo = AT + X{z)AT^ H , with X{z) = g^osiz) /2{1 + z)^, where AT = {1 + z)Ate is the arrival time neglecting 

the cosmic acceleration/deceleration (the second term). Now for AT ^ 10^ sec, the time lag of the arrival time due 
to the cosmic acceleration/deceleration amounts to the order of second ~ 10^^/(3 x 10^^) ~ 1 [sec]. From Eq. (|^), if 
X{z) is positive, then dta{to) > dta{te)- This clearly means that our universe is accelerating. Therefore the value of 
this time lag is the direct evidence for the acceleration/deceleration of the universe. 

As shown in if the accurate value of X{z) at a single point Zs is available it is possible to determine Vlka as 
= {1 ~ {dr{zs)/dz)'^{l + Zs)'^{Hq — 2X{zs))^}{r{zs)^ Hq}~^ , where we have assumed that the quantity r{z) = 
dL{z)/{l + z) is obtained accurately, e.g., by SNAP. Even if the accurate values of X{z) are not available for any 
points, we may apply the maximal likelihood method to determine ^ko- Using the value of O^q thus determined, we 
can obtain the equation of state of our universe without any theoretical assumption on its matter content ||^ . 

B. Evolution of Phase of Gravitational Waves from Coalescing Binary at Cosmological Distance 

Let us study an inspiraling compact binary system that evolves secularly by radiating gravitational wave 
For simplicity we study a circular orbit and evaluate the gravitational wave amplitude and the energy loss rate by 
Newtonian quadrupole formula. We basically follow analysis of Cutler & Flanagan but properly take into account 
of effects of accelerating motion. The Fourier transform h{f) — e^^^^'^h{t)dt for the wave h{t) is evaluated using 

the stationary phase approximation as h{f) = KdL{z)~^Mc^^ f^'^^^ exp[i$(/)], where K is determined by the angular 
position and the orientation of the binary relative to the detector, and Ale is the chirp mass of the system. Keeping 
the first order term of the coefficient X{z), the phase <&(/) of the gravitational wave becomes 
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$(/) = 2nfU -^c-^ + -(87rM,./)-5/3 - 



X(z)/-i3/3Af-io/3^-i3/3, 



(2) 



where tc and 0c are integral constants and Mcz — Mc{l + z) is the redshifted chirp mass 
If we include the post-Newtonian (PN) effects up to P^'^N-order, the term 3/4(87r Mcz/) 
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modified as USnMczfy^/^ 
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in eq.(y) should be 

) x + (4/3 - 167r)x3/2 + • • • , where x = (7rMJ(I + z))2/3 = 0{v^/c^) 



9 ^336 ' 4Mc 

is the PN expansion parameter with Mt being the total mass of the binaries. The term proportional to /3 in P^-^N 
order (oc a;^/^) is caused by the spin effect ||l]] , |l^ . In general P^N contribution depends on the frequency / as 
0(/(^^+^^^/3) and is largely different from the dependence /"^^/'^ caused by the cosmic acceleration. This difference 
is very preferable for the actual signal analysis. 



C. the estimation error 

For the circular orbit of the binary neutron stars (NSs) of mass Mi and M2 with the separation a at the redshift 
z, the frequency of the gravitational waves / is given by / = O.IHz(l + z)^^(Aft/2.8MQ)^/^(a/15500km)~3/2^ The 
coalescing time tc, the number of cycles N cycle and the characteristic amplitude of the waves he are given by 

tc = 7(1 + z)(Afi/1.4M0)~i(M2/1.4M0)-i(Mt/2.8MQ)-i(a/155OOkm) V (3) 
Ncycie = 1.66 X lO^(Mi/1.4M0)^i(M2/1.4M0)-i(Mt/2.8M0)-i/2(a/155OOkm)5/2 (4) 
he = 1.45 X 10-^^(1 + z)5/^(Mc/1.2M0)^/S(//O.lHz)-i/^(dL/lOGpc)-i (5) 

Let us evaluate how accurately we can fit the parameter X{z). We take six parameters Xi — 
{A, Mcz, fJ-zjtc, (f>c, Mcz^'^^^X{z)} in the matched filtering analysis up to IPN-order for the phase ^{f) and New- 
tonian order for the amplitude Here A is the amplitude of signal Kd^^z)^^ Mc^^ in the previous subsection 
and fiz is the redshifted reduced mass fiz — + z)MiM2/Mt. As the chirp mass Mcz can be determined quite 

accurately, we simply put AX{z) — A ^Mcz^^^^ X{z)^ jMcz^^'^ . For simphcity we fix the redshift of sources at z = 1 

and calculate S/N and the error AX for equal mass binaries with various integration time At before coalescence. We 
use the effective factor 1 / \/5 for reduction of antenna sensitivity due to its rotation ||l[ . For the present analysis we 
neglect the binary confusion noise since double White Dwarf binaries do not exist at frequency / ^ O.lHz fl^ . 

We found that we can detect NS-NS binaries at z = 1 with S/N ~ 20000 and AX/t^^ ~ 7.0 x 10~3 for integration 
time At = 16yr {Ncyde ^ 10^ orbital cycles), and S/N ~ 10000 and AX/tf^^ ~ 1.26 for At = lyr. With this detector 
it would be possible to determine X{z) and obtain the information of the cosmic acceleration quite accurately. With 
AT = 16yr we have the estimation error for the redshifted masses as AMcz/Mcz = 1-5 x 10^^^, A/i^Z/Zz = 4.2 x 10~® 
and for the wave amplitude AA/A ^ (S/N)^^ — 5x 10^^. Although the more detailed study is needed to estimate the 
error of the binary inclination angle, it is expected that the luminosity distance can be determined accurately so 
that the redshift z can be determined using the inverse function z — (distance) of the accurate luminosity distance 
from e.g. SNAP. As a result we can know two (not redshifted) masses Mi and M2 for ~ lO'^ binaries per year up to 
z = 1 . This will be large enough to establish the mass function of NS which would bring us important implications 
for the equation of the state of the high density matter and the explosion mechanisms of Typell supernovae. 

As the S/N and the estimation error scale as S/N oc h^,}^^ and AX oc hrms, we can attain AX/tQ^ ~ 7.0 for the 
integration time T — 16yr using a less sensitive detector with hrms ~ 10~^^ (1000 times worse). Even though the 
error bar AX is fairly large for this detector, the likelihood analysis would be an efhcient approach to study the cosmic 
acceleration. Considering the estimated cosmological coalescence rate of NS-NS binaries (<] 2 x (10Gpc/350Mpc)3 ~ 
lO^yr^^) |lj], we may expect the decrease of the estimation error AX roughly by a factor of ^ 1/300 ~ 1/VTo^. 



D. Acceleration in the Very Early Universe 

In the inflationary phase there was an extremely rapid acceleration of the universe. In this phase the gravitational 
waves were generated by quantum fluctuation ||l5|. With CMB quadrupole anisotropics measured by COBE, the 
slow-roll inflation model predicts a constraint on the stochastic background flaw ^ 10^^^ — 10~^^ at / '--^ O.lHz [ p^ . 
Ungarelli and Vecchio |]l^ discussed that the strain sensitivity hrms 10"^** is the required level at f ^ O.lHz for 
detecting ^Igw 10~^^ by correlating two detectors for decades (see also Ref. jl^). It is important to note that the 
band / > 0.1 Hz is free from stochastic backgrounds generated by White Dwarf binaries. The radiation from neutron 
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stars binaries is present in this band and it is indispensable to remove their contributions accurately from data stream, 
where effects of the cosmic acceleration would be highly important. Thus measurement of the present-day cosmic 
acceleration is closely related to detection of the primordial gravitational wave background that is one of the most 
interesting targets in cosmology. If DECIGO with hrms ^ 2 x 10~^^ at / ~ O.lHz is available we can detect the 
primordial gravitational waves background even if the energy density is extremely low flaw ^ 10"^'' by correlating 
two detectors for a decade. 

Confusion noise due to NS-NS (or NS-BH, BH-BH) binaries might be important in the band / ^ O.lHz. Ungarelli 
and Vecchio [ pT[ investigated the critical frequency fg where we can, in principle, remove signal from individual NS- 
NS binaries by matched filtering analysis and the observed window becomes transparent to the primordial stochastic 
background. They roughly estimated fg ~ O.lHz where the number of binaries per frequency bin (~ 10~®Hz) is less 
than one. But binaries around / '^-^ /g ~ 0.1 chirp significantly within observing time scale and the situation would 
be more complicated than monochromatic sources ||l^. Although more detailed analysis is needed, a much smaller 
NS-NS coalescence rate than ~ lO^yr"-'^ might be required for our analysis to be valid. 



IV. DISCUSSIONS 

The determination of the angular position of the source is crucial for matching the phase [Q . The phase modulation 
at the orbital radius lAU corresponds to 2AU/c ~ 1000[sec]. Thus, in order to match the phase within the accuracy 
of 0.1[sec] we need to determine the angular position with precision ~ 0.1/1000 [rad] ~ 20". In the matched filtering 
analysis we can simultaneously fit parameters of the angular position as well as the relative acceleration between the 
source and the barycenter of the solar system. Due to their correlation in the Fisher matrix, the measured acceleration 
would be somewhat degraded if we cannot determine the angular position by other observational methods. Using 
the gravitational wave alone, we can, in advance, specify the coalescence time and the angular position of the source 
within some error box. If coalescence of NS-NS binaries would release the optical signal {e.g. Gamma Ray Bursts as 
proposed by we may measure the angular position accurately by pointing telescopes toward the error box at the 
expected coalescence time from the chirp signal. Therefore we have not tried to fit the angular position of the source 
in the matched filtering method ||l[ . We might also determine the redshift of the source by using optical information 
of host galaxies. 

Let us discuss the effects of the local motion giocai of the emitteron the second derivative cPto/dt^ (see e.g. Ref. po|). 
As the effect of bulk motion of galaxy is much smaller than cosmological effect, we estimate the internal acceleration 
within the galaxy based on the observational result of NS-NS binary PSR 1913-1-16. As shown in Table 1 of [po[ , 
the dominant contribution of its acceleration x comes from the global Galactic potential field and has time scale 
c/x ~ 10io(^e/10kpc)(T/rot/200kms"^) that can be comparable to the cosmic signal gcos where Re is the effective 
radius of the acceleration and Vrot is galactic rotation velocity. However the contamination of local effect giocai can 
be reduced by taking the statistical average of many binaries as {gcos giocai) — {gcos)- We also note that the 
cosmological change in phase of a coalescing binary [given by the last term in Eq. (3)] may have other applications, 
and may under certain circumstances be observable by the planned LISA mission. 

In conclusion we would like to encourage the further design study of DECIGO and the theoretical study of the 
sources of gravitational waves for DECIGO. Even if we may not see the construction of DECIGO in our life since the 
highly advanced technology is needed, we are sure that our children or grandchildren will decide and go DECIGO. 
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would like to thank H. Tagoshi for useful discussions and M. Sasaki for his invaluable advice. This work was supported 
in part by Grant-in- Aid of Monbu-Kagaku-Shyo Nos.11640274, 09NP0801 and 0001416. 
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FIG. 1. Sensitivity (effectively S/N=l) for various detectors (LISA, DECIGO, LIGOII and a detector 10^ times less sensitive 
than DECIGO) in the form of hrms (solid lines). The dashed line represents evolution of the characteristic amplitude he for 
NS-NS binary at z = 1 (filled triangles; wave frequencies at lyr and 10 yr before coalescence). The dotted lines represent the 
required sensitivity for detecting stochastic background with flaw = 10~^® and Qaw = 10"^" by ten years correlation analysis 
(S/N=l). 
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